Composite Pairings in Chirally Stabilized Critical Fluids 
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We study a one-dimensional electron gas in a special antiferromagnetic environment made by two 
spin- 1/2 Heisenberg chains and the one-dimensional two-channel Kondo-Heisenberg lattice away 
from half-filling. These models flow to an intermediate fixed point which belongs to the universality 
class of chirally stabilized liquids. Using a Toulouse point approach, the universal properties of the 
models are determined as well as the identification of the leading instabilities. It is shown that these 
models exhibit a non- Fermi liquid behavior with strong enhanced composite pairing correlations. 
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The concept of odd-frequency superconducting order- 
ing was conceived by Berezinskii ]lj who considered spin- 
triplet odd-frequency pairing as an alternate of the con- 
ventional theory 3 He pairing. The possibility of such 
unusual pairing was revived nearly one decade ago by 
Balatsky and Abrahams Q in the context of singlet su- 
perconductors and by Emery and Kivelson |3j in the 
two-channel Kondo problem with the emphasis laid on 
the composite nature of the order parameter for odd- 
frequency pairing. This idea of composite pairing was 
further developped by Coleman et al. [Q who proposed 
that heavy-fermion superconductors might involve odd- 
frequency triplet pairing. 

One of the most important difficulty in finding some 
specific lattice models for realization of odd-frequency 
pairing stems from the fact that it requires a controlled 
analysis in the strong coupling regime. Such an ap- 
proach is possible in some extreme limits as in one di- 
mension where non-perturbative techniques are available. 
Strong odd-frequency singlet pair correlations were, in 
particular, identified within the bosonization approach 
H of the one-dimensional single channel Kondo lattice 
in an anisotropic limit (^J^). Another candidate is the 
one dimensional single channel Kondo-Heisenberg lattice 
(KHL) which consists of a one dimensional electron gas 
(1DEG) interacting with an antiferromagnetic Heisen- 
berg spin- 1/2 chain by a Kondo coupling. Away from 
half filling, this model has a spin gap || and is thus ex- 
pected to exhibit enhanced pairing correlations. In par- 
ticular, it has been shown recently Q that the dominant 
instabilities are of a composite nature in a certain regime 
of the parameters. In this letter, we shall study two dif- 
ferent generalization of the one dimensional KHL with 
two channels which have no spin gap but exhibit domi- 
nant unconventional pairing instabilities. A first gener- 
alization is to consider a 1DEG coupled by a Kondo cou- 
pling to two non-interacting antiferromagnetic Heisen- 
berg spin-1/2 chains. This model belongs to the more 
general class of Luttinger liquids in active environments 
introduced in the context of the striped physics [jio] |l2| . 
A second generalization consists of conduction electrons 



with a two-fold orbital degeneracy interacting with a pe- 
riodic array of localized spins i.e. the two-channel KHL. 
It has been shown in Ref. jl3| that these two models 
away from half filling exhibit a critical phase governed 
by a fixed point which belongs to the class of chirally 
stabilized fluids In this work, we shall compute 

correlation functions of physical observables by means 
of a Toulouse point solution and characterize the dom- 
inant instabilities of both models which turn out to be 
of a composite nature. Very recently, the one dimen- 
sional N-channel KHL has been investigated by Andrei 
and Orignac [jl5| using a conformal field theory (CFT) 
approach. These authors found that the leading insta- 
bility is of a composite pairing type when N < 4 Jl5| . 
In this respect, the exact solution at the Toulouse point 
provides an independent and physically transparent ap- 
proach for the particular 2-channel KHL. 

The models. The first model (referred in the following 
as a 1DEG in a special active environment) consists of 
a 1DEG away from half filling coupled symmetrically by 
a Kondo coupling (Jk > 0) with two non-interacting 
antiferromagnetic (Jh > 0) spin-1/2 Heisenberg chains: 
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Here, Ci a is the conduction electron operator at site i with 
spin a =t, L S C i — c\ a a a fjCifjj2 (a being the Pauli matri- 
ces) stands for the electron spin operator at site i whereas 
the localized spin operator at site i on the chain of index 
a = 1,2 is denoted by S a i- In the limit Jk *C £, Jff, 
the low energy behavior of the model can be determined 
using the continuum limit of the electron operator in 
terms of right and left moving fermion fields: Ci a — » 
R a e lkFX +L a e~ zkFX , x = iao (hp being the Fermi momen- 
tum and ao the lattice spacing). The continuum descrip- 
tion of the conduction spin density operator S c can then 
be expressed in terms of a spin current 3qr,l belonging to 
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the SU(2)i Kac- Moody (KM) algebra [§ and a bosonic 
field <& c = &cL + &cR that accounts for the charge degrees 
of freedom: S c — > 3 0R + Jql + cos{2kpx + ^/2tt<& c )yi , n 
being the staggered magnetization. In the same way, the 
spin densities of the surface chains are represented as 
(see Refs. 05)): S a (x) = 3 a (x) + (-l) x /°°n a (x) where 
Ja = 3 aR + JaL and n a are respectively the uniform and 
staggered parts of the magnetization. The chiral spin 
currents 3 a R,L belong also to the SU(2)i KM algebra. 
With this low energy description, the Hamiltonian den- 
sity of the lattice model (jl]) for incommensurate filling 
reads as follows in the continuum limit: 
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where we have neglected all oscillatory, marginal irrele- 
vant contributions as well as current-current interactions 
of the same chirality. As shown in Ref . |D| , at the strong 
coupling fixed point, these latter interactions only lead to 
rcnormalization of velocities and logarithmic corrections. 
In Eq. (H), 9 C = — & cR is the dual charge field, vf 
the Fermi velocity, vh the spin velocity of the magnetic 
environment, and g ~ Jk&o > is the interacting cou- 
pling constant; l R . L = 3 1RtL + 3 2 r,l is a SU(2) 2 KM 
current being the sum of two SU(2)i KM currents. 

The second model considered in this letter is the one- 
dimensional two-channel KHL: 
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where now the electron operator Ci aa has a channel index 
a = 1,2 and interacts by a Kondo coupling {Jk > 0) 
with a periodic array of localized spins So^. Away from 
half filling, the continuum limit of the Hamiltonian (|^) 
proceeds in the same way as in the first model and with 
the same degree of approximations than in Eq. (|J) , one 
has: 
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where $ ac is a bosonic field (9 QC being the dual field) 
associated with charge fluctuations in the ath channel. 
The chiral SU(2)i currents J a fl,i correspond to the chiral 
uniform part of the electron spin density whereas Jo_r,l 
are the chiral SU(2)i currents of the local moments. 

Toulouse point solution. The next step of the ap- 
proach is to use the representation of two SU(2)i cur- 
rents in terms of four Majorana fermions £° and £ [p|,p7|: 



Ir,l — — i A £/?„l/2, Jir,z — 3 2R .l — ^r,l£, r , l - 
The Hamiltonians (@,||) share then the same structure: 
Hi, 2 = H° c + H°(f) + H where H° c is the free Hamil- 
tonian of the charge bosonic fields and 7i°(£°) is the 
free Hamiltonian of the massless Majorana fermion £°. 
All nontrivial physics is incorporated in TL which sepa- 
rates into two commuting and chirally asymmetric parts: 
H = Hi + H2, [Wi,W a ] = 0with 
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with v\ = vh,vq — vf for the 1DEG in a special active 
environment and the reverse for the two-channel KHL; 
TL 2 is obtained from TLi by inverting chiralities of all 
the currents. For an antiferromagnetic Kondo coupling 
(<? > 0), the interacting part of TL is marginally rele- 
vant. Naively one will thus expect that the model will 
enter a massive strong coupling region. However, due 
to the chiral asymmetry of TL, the effective interaction 
flows towards a conformally invariant intermediate fixed 
point with a smaller central charge than in the ultravi- 
olet. Stated differently, the interaction produces a mass 
gap in some sectors but there are still some degrees of 
freedom that remain critical in the infrared (IR). The 
nature of the fixed point of chirally asymmetric current- 
current models with the structure (^|) have been deter- 
mined by Andrei et al. by a combination of CFT and 
Bethe ansatz techniques. This fixed point belongs to the 
class of chirally stabilized fluids introduced in Ref. 
describing the IR behavior of one-dimensional interact- 
ing chiral fermions. A simpler way to identify the critical 
degrees of freedom of model (||) in the IR limit can be 
done by means of a Toulouse point approach as in the 
two-channel Kondo model ||. By allowing anisotropic 
interaction (g — > <?j|,<7j_), the U(l) version of the Hamil- 
tonian TL can be mapped onto free Majorana fermions for 
a special value of the interaction gjj . Although the posi- 
tion of the solvable point is non-universal, the Toulouse 
limit solution captures the physical and universal proper- 
ties of the models including the SU(2) case. The details 
of the approach can be found in Refs. [ [T8||l5| ] and we shall 
now briefly review it to fix the notations. The starting 
point of the Toulouse solution is the Abelian bosoniza- 
tion of all SU(2)i currents 3 a ,a — 0,1,2 in terms of 
massless bosonic fields ip a ||: J^ RL = {l/V2~Tr)d x ip a R,L, 
J+ RL = (l/2-n-a )e^ i ^t'> R ' L . Introducing the sym- 
metric combination of the fields: ip + — (<pi + tp2)/y2, 
the SU(2) 2 current I writes: I RL = (l/y/TF)d x <p+ R ,L, 

I R L = (i/y/TTOo)^ R L Ke Tl ^^ ip + Ii - L . A fermionic zero- 
mode operator n has been introduced to ensure the 
proper anticommutation relations. The following canon- 
ical transformation is then performed: 
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where $o and (respectively 6i and 62) are the dual 
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fields associated with tpo and y>+ (respectively $1 and 
i> 2 ). For a special positive value (Toulouse point) of 
<7||, = 47r(w + vi)/3), the arguments of the inter- 
acting terms become those of free fermions so that they 
can be refermionizcd further with the introduction of a 
pair of Majorana fields (77, £) using the correspondence: 
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where to = g±/2irao, u± = (2v± — vq)/3, and u 2 = 
(2v Q — ui)/3. The Toulouse solution is stable provided 
that all velocities are positive i.e. 1/2 < vo/t'i < 2. 
The first two terms in Eq. (|t]) describe decoupled free 
massless bosonic and Majorana fields, <§i and £, con- 
tributing to criticality with the central charge: c = 3/2. 
The remaining part of the Hamiltonian (^) has a spec- 
tral gap to and describes hybridization of the Majorana 
£ 3 and r] fermions with different chiralities. Adding the 
contribution of the excitations that do not participate 
in the interaction, the total central charge in the IR of 
the 1DEG in a special active environment (respectively 
the two-channel KHL) is c = 3 (respectively c = 4). 
Apart from the charge degrees of freedom, the elemen- 
tary excitations of the models in the IR are of different 
nature: The magnetic excitations are spinons defined as 
\Jir/2 kinks of the bosonic field $1 describing an effective 
S=l/2 Heisenberg chain and nonmagnetic singlet excita- 
tions associated with the two massless Majorana fermions 
£°, C (critical Ising degrees of freedom) referred as pseu- 
docharge degrees of freedom in Rcfs. [ pj||l3| ]. These Z 2 
excitations play a crucial role in the nontrivial IR physi- 
cal properties of the models as we shall see now. 

Physical properties of the 1DEG in a special active en- 
vironment. The Green's function for the right-left mov- 
ing fermions can be computed using the bosonic represen- 
tation: (R, L) a = (K CT /V27rao) e lT ^/ i e ±{ ^^'' R - L where 
r T ,j. = ±1, ip aR . L = ($ c r,l + T a ip aRtL )/V2, and k ct are 
Klein-factors to insure the correct anticommutation rela- 
tions between fermion fields of different spin index [fl9|| . 
Using the Toulouse basis (||), one then obtains at the IR 
fixed point the estimate: 
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which means that the system displays non-Fermi liquid 
properties; the left-moving Green's function is obtained 



from Eq. ( 
teresting p 



) by the substitution: i — ► — i. The most in- 
lysical quantities are the correlation functions 



correlation functions have been computed in Ref. |T3|] and 
the slowest ones are the staggered parts that decay as 
^,-3/2 u s i n g the previous Abelian bosonization of the 
chiral fermions and the Toulouse basis ([|), we find in 
the electronic sector the following representation for the 
charge density wave (CDW), spin density wave (SDW), 
singlet (SS) and triplet (ST) superconducting order pa- 
rameters in terms of the different critical fields at the IR 
fixed point: 
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Oss = -iL a (* v ) a0 R ~ e"^ ' cos (v^i) CaCl 
Ots = -iL a (<?<r») a/3 Rp * ie- iv ^ 6c Cfla 
cos rV27r6i) , - sin (V^tQi) , sin (v^ii) • (9) 

All these conventional order parameters have thus the 
same scaling dimension (2) and the corresponding pair 
correlation functions decay as x~ A i.e. much faster 
than in the one dimensional metals. One should notice 
that the Z2 pseudocharge excitations contribute in (||) 
through the density energy operator i(rCl- In a critical 
Ising model, there are other primary fields (order and dis- 
order parameters) that have a smaller scaling dimension 
(1/8) which are highly nonlocal in terms of the original 
Majorana fermion M. This leads us to investigate the 
possibility that the ground state might be characterized 
by composite order parameters: staggered odd-frequency 
singlet pairing (c-SP) and staggered composite CDW (c- 
CDW). At the IR fixed point, we find the following cor- 
respondence for these order parameters: 
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where /^0jM4 (respectively (Jq,<J4) are the Ising disorder 
(respectively order) parameters associated with the Ma- 
jorana fermions £°,C- /.From Eq. (fiOj), we deduce the 
leading asymptotics of the correlation function corre- 
sponding to the composite order parameters: 
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of the various possible order parameters. The spin-spin 



so that the composite order parameters induce the domi- 
nant instabilities and have enhanced long-range coher- 
ence. These fluctuations can be made even more en- 
hanced upon switching on a very small ferromagnetic 
(.912 < 0) exchange interaction between the spin chains: 
012 = 512S1 • S 2 . Indeed, at the IR fixed point, this per- 
turbation mostly affects the pseudocharge sector leaving 
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intact the magnetic one: O12 — 9i2^r^l- The Majorana 
fermion £° acquires a positive mass i.e. its associated 
Ising model is in the disorder phase: (fj,o) ^ and the 
Z2 symmetry with respect to the interchange of the two 
chains is now broken. All other degrees of freedom remain 
critical so that the pair correlation of the conventional or- 
der parameters (^|) still decay as a;" 4 whereas from Eq. 
( |i"o| ) one immediately observes that the composite ones 
decay now as x~ 5 / 4 instead of x~ 3 / 2 at 912 = 0. 

The two-channel KHL case. The determination of the 
ground-state physical properties of the two-channel KHL 
proceeds in the same way as in the previous model. We 
first use the Abelian bosonisation of the right-left moving 
fermions: [R,L) aa = {n aa /^/2^) e ^/i e ±i^ VaaR , L 
where ip a aR,L = {^acR.L + Ta^aH.h)! ' V%, and 

^aa are 

some Klein-factors jF9|. The Green's functions of the 
chiral fermions have a similar structure than in Eq. 
(^) apart from they fall now with the distance with 
a power 5/4. The leading asymptotics of the corre- 
lation function of the localized spins have been com- 
puted in Ref. [T^ ] and in particular the staggered part 
decays as x~ 3 . Using the Toulouse basis (|g) and in- 
troducing the charge and relative charge bosonic fields 
®±cR,l = {®icR,L ± $2c.r,l)/\/2, the order parameters 
for the electronic degrees of freedom can also be expressed 
in terms of the critical fields at the IR fixed point. We 
find the following correspondence for the conventional or- 
der parameters: 
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where 1$ = /io/M cos (y/ir$- c ) + ia^aA sin (^/tt^-c) and 
I® is obtained from Z$ by replacing <I>_ C by its dual 
field 0_ c . The scaling dimension of all these conven- 
tional order parameters is 5/4 and their corresponding 
pair correlation functions decay thus as x~ 5 ^ 2 i.e. have 
no enhanced long-range coherence. Notice, however, that 
this decay is much slower than in the 1DEG in a special 
active environment case. Similarly, we find the follow- 
ing description at the IR fixed point for the staggered 
composite order parameters: 
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so that the scaling dimension of these operators is 3/4 
and the pair correlation function of the composite order 



parameters is still given by Eq. (0). The dominant 
instabilities of the two-channel KHL are thus of a com- 
posite nature as in the 1DEG in a special active environ- 
ment. Finally, one should note that the scaling dimen- 
sions of the physical observables of the two-channel KHL 
obtained by the Toulouse point solution coincide with 
those of the recent CFT approach In this respect, 
we stress that, in contrast with the latter approach, the 
Toulouse point solution unables us to take into account 
the inherent velocities anisotropy in the problem. 

In summary, we have shown that the two generaliza- 
tions of the one dimensional KHL considered in this let- 
ter exhibit a nontrivial non- Fermi liquid low-temperature 
phase belonging to the class of chirally stabilized fluids 
with strong enhanced staggered composite pairing corre- 
lations. 
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